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Spectral statistics of molecular resonances in erbium isotopes: How chaotic are they? 
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We perform a comprehensive analysis of the spectral statistics of the molecular resonances in 
^®®Er and ^®®Er observed in recent ultracold collision experiments [Erisch et ah, Nature 507 , 475 
(2014)] with the aim of determining the chaoticity of this system. We calculate different independent 
statistical properties to check their degree of agreement with random matrix theory (RMT), and 
analyze if they are consistent with the possibility of having missing resonances. The analysis of the 
short-range fluctuations as a function of the magnetic field points to a steady increase of chaoticity 
until B ~ 30 G. The repulsion parameter decreases for higher magnetic fields, an effect that can 
be interpreted as due to missing resonances. The analysis of long-range fluctuations allows us 
to be more quantitative and estimate a 20 — 25% fraction of missing levels. Finally, a study of 
the distribution of resonance widths provides additional evidence supporting missing resonances of 
small width compared with the experimental magnetic field resolution. We conclude that further 
measurements with increased resolution will be necessary to give a final answer to the problem of 
missing resonances and the agreement with RMT. 
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I. INTRODUCTION 

Understanding the structure of a molecule is a complex 
task that usually requires input from both experimental 
spectral data and electronic structure calculations. This 
joint effort has the goal of identifying the relevant sym¬ 
metries and quantum numbers, ultimately allowing one 
to translate an initial apparently random collection of 
spectral lines into a coherent set of vibrational bands, 
rotational lines, etc BE], and is analog to the process 
necessary to understand atomic or nuclear spectra. 

In atomic and, especially, nuclear spectroscopy, it has 
long been recognized that some systems feature very com¬ 
plex spectra which defy this procedure. In essence, inter¬ 
actions between the different particles in the system are 
so strong that the many-body system becomes strongly 
correlated. The analysis of its spectrum can not iden¬ 
tify each line unambiguously with a dominant set of 
quantum numbers, and one resorts instead to statisti¬ 
cal studies. Such chaotic quantum systems present very 
specific statistical properties in their spectra and wave 
functions that can be traced to results from random ma¬ 
trix theory (RMT) [3J|3]. A limiting case of particular 
relevance holds when the quantum system has a fully 
chaotic classical analog. In this case, according to the 
original Bohigas-Giannoni-Schmit (BGS) conjecture, its 
spectrum will show “the same fluctuation properties as 
predicted by the Gaussian orthogonal ensemble (GOE)” 
[5|. This conjecture has now been given firm grounds at 
the semiclassical level [5]. 

There is an enormous amount of numerical and ex- 
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perimental evidence in support of the relationship be¬ 
tween the properties of chaotic systems and the prop¬ 
erties of random matrices, both in systems with classi¬ 
cal analog and in many-body quantum systems la [Tj. 
The most stringent tests so far have been performed on 
large data sets of nuclear levels Big as well as on mi¬ 
crowave billiards [3]. There exist also several studies on 
the quantum-chaotic character of small molecules spectra 
from stimulated-emission pumping spectroscopy mm- 
In the last decade, atomic and molecular physics has been 
enriched by new experimental and theoretical tools for 
high-resolution spectroscopy at very low energies mm 
as well as novel control and measurement protocols |13T 
m- In this context, the field of quantum chaos has re¬ 
cently received new input from ultracold collision exper¬ 
iments m- In this work we present a detailed study of 
these new data. 

Analyzing experimental systems to check whether they 
behave quantum chaotically or not is far from a straight¬ 
forward process as several effects may lead to deviations 
from RMT predictions. First, many systems present in¬ 
termediate behavior between regularity and chaos. Sec¬ 
ond, the statistical analysis to compare with RMT re¬ 
sults should be done for complete sequences of levels with 
the same quantum numbers. There have been a num¬ 
ber of papers on how to analyze incomplete or imper¬ 
fect spectra [ini [TM^T)] . They study both short and long 
range spectral fluctuations. Here, by incomplete spectra 
we mean that there is a fraction of missing levels unac¬ 
counted for experimentally, and by imperfect spectra we 
mean that there are wrong assignments of quantum num¬ 
bers in the experimental spectra. Crucially, the effects of 
both missing levels and imperfect spectra are similar on 
short-range fluctuations, reducing the level repulsion and 
making them appear non-chaotic. On the other hand, 
their effects are very different for long range fluctuations, 
which can be then used to discern these two effects. For 
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example, Ref. |20] showed how the power spectrum is 
a useful statistic to discriminate between them, noting 
that missing levels affect low frequencies (i.e., correla¬ 
tions with longer range) first, while mixed symmetries 
affect high frequencies (correlations with shorter range) 
first. 

Very general formulas can be derived for the power 
spectrum as a function of the number of mixed symme¬ 
tries, the number of levels in each symmetry, and the 
fraction of missing levels for each symmetry. In practice, 
in order to have simple formulas useful to extract the 
fraction of missing levels or the number of mixed symme¬ 
tries from the experimental data, assumptions about the 
original spectra are needed and the mentioned references 
assume that the original spectra follow the predictions 
of RMT, i.e., that the underlying system is fully chaotic. 
What would happen if this were not the case and the 
original spectra had intermediate statistics? Likely, one 
will obtain a result closer to Poisson statistics (i.e., reg¬ 
ular dynamics) than what the chaoticity of the original 
spectra would suggest and that will be very difficult to 
interpret. However, the fact that short- and long-range 
correlations give mainly independent information can be 
used to our advantage. 

Arguments for the study of different statistical proper¬ 
ties comparing short- and long-range level statistics with 
wave function statistics have been given before, for exam¬ 
ple noting that the transition between regular and chaotic 
spectra is not universal and can reflect in different ways 
in the different statistics [5T]. Incomplete and imper¬ 
fect spectra are other instances when complementary ap¬ 
proaches are needed. For example, wave-function statis¬ 
tics give another completely independent test of the data 
versus RMT predictions and, if the missing levels are the 
ones with smaller widths, as it is expected for finite exper¬ 
imental resolution, one can obtain this fraction from the 
analysis of the statistical distribution of level widths in 
comparison with the Porter-Thomas prediction of RMT 

m- 

The general strategy to follow is, then, to study as 
many independent statistical properties as possible and 
not restrict oneself to only short-range (nor only long- 
range) correlations. If the complete analysis gives a con¬ 
sistent value for the number of missing levels, one can 
conclude that this number reflects the experimental re¬ 
sults and that the original data do follow RMT. Other¬ 
wise, several effects may be playing a relevant role, and 
the original data probably follow intermediate statistics. 

Here, we apply such a comprehensive approach to the 
data on molecular resonances in a lanthanide atom re¬ 
cently reported by Frisch et al. ffH- The authors of that 
work studied collisions of trapped cold erbium atoms as 
a function of the magnetic field m- The Feshbach reso¬ 
nances of bosonic isotopes ^®®Er and ^®®Er were studied 
with high-resolution trap-loss spectroscopy. A statistical 
analysis was made of the positions of the resonances using 
the tools of RMT. The behavior obtained was intermedi¬ 
ate between Poisson and Wigner-Dyson (WD) statistics 


although somewhat closer to WD. From these results it 
was concluded that the dynamics of collisions and for¬ 
mation of Er 2 molecules is very complex and presents 
characteristics of quantum chaos. We have extended the 
analysis of the experimental data in Ref. perform¬ 
ing a study of short and long-range fluctuations in level 
positions as well as of level widths. 

The paper is organized as follows. In Sec. [H] we present 
a detailed analysis of the short and long range fluctua¬ 
tions of level statistics. Our conclusions here agree gener¬ 
ally with those in Ref. m- At the same time, our results 
enable an analysis of the magnetic field dependence of the 
degree of chaoticity, leading us to a quantitative account 
of the origin of the deviations from Wigner-Dyson of both 
the short- and long-range correlations in the spectra. In 
addition, our calculations allow to estimate a fraction of 
missing levels at high magnetic fields. Then, in Sec. [m] 
we present calculations of the level width statistics. From 
this study, it follows that there is a substantial fraction of 
missing levels with small widths—i.e., with widths below 
the current magnetic field resolution—which precludes 
reaching definitive conclusions on the chaoticity of the 
system. The conclusions from these studies are summa¬ 
rized in Sec. CYl 


II. STATISTICS OF ENERGY LEVELS 
A. Unfolding the spectra 

The unfolding of the spectra is the procedure by which 
the eigenenergies of the system are transformed to 
new unfolded energies so the density of levels is lo¬ 
cally constant and equal to unity [3]. It is a heuristic 
procedure, essential in order to compare the statistical 
analysis of the spectra with RMT results. Doing the un¬ 
folding erroneously can give rise to misleading signatures 
of chaos [23]. Fortunately, for the erbium isotopes it is 
relatively uncomplicated. In Ref. m it was noticed that 
the density of resonances was constant for magnetic fields 
B > 30 G. The unfolding in that case is trivial and the 
results for the spacing distribution and other statistical 
measures very reliable. 

We have chosen to use two unfolding methods to cross¬ 
check the robustness of our results. The first one is a local 
unfolding where an average density of states is calculated 
in a window of energies (in the present case, a window 
of magnetic fields) and this window is moved as we un¬ 
fold the spectra |3]. This local unfolding gives good re¬ 
sults for the short-range correlations such as the nearest- 
neighbor spacing (NNS) distribution, but it underesti¬ 
mates dramatically long-range correlations [23] . The sec¬ 
ond method is a polynomial fit of the accumulated den¬ 
sity of states N{B )—also referred to as the “staircase 
function”—as a function of external magnetic field, which 
is a measure of the resonance excitation energy. The sim¬ 
plest case of a constant density of states amounts to a 
linear fit to the accumulated density of states. In Fig. 
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FIG. 1. (Color online) Comparison of the accumulated num¬ 
ber of levels as a function of the magnetic field N{B) for ^®®Er. 
Symbols represent data from Ref. m while lines are linear 
(dashed) and polynomial (full) fits. Similar results hold for 
^®®Er (not shown). 


we show a comparison of the actual accumulated density 
of states with the smooth approximation to a fifth order 
polynomial for all data 0 < B < 70G, and a linear fit to 
data for B > 30G for both isotopes. The results of the 
polynomial fit are very good for all magnetic fields. 


B. Short range correlations 

Once a reliable unfolding has been made we can start 
studying spectral statistics. The most used tool for com¬ 
paring with RMT results is the nearest-neighbor spacing 
(NNS) distribution, P{s), that measures short range cor¬ 
relations and, in particular, the level repulsion between 
spectra. The NNSs are defined as the differences between 
consecutive unfolded energies Si = e^+i — Ci, so that the 
mean (s) = 1. We have checked that both local and 
polynomial unfoldings give essentially the same results 
for P{s). 

In order to study level repulsion we fit the results of the 
P{s) to the Brody distribution [24], as was made in the 
original experimental paper The Brody distribution 
is a phenomenological distribution that interpolates as a 



FIG. 2. (Color online) Brody parameter vs. magnetic field 
from fits of the nearest-neighbor spacing distribution P{s) on 
a moving window of 60 levels for ^®®Er (empty red circles) 
and ^®*Er (filled blue squares). The dashed horizontal lines 
indicate the limiting values corresponding to Poisson and WD 
distributions. The vertical error bars indicate the errors in the 
fit. 

function of the repulsion parameter f3 between the Pois¬ 
son distribution corresponding to regular spectra (when 
/? = 0) and the Wigner-Dyson (WD) distribution corre¬ 
sponding to chaotic spectra {j3 = 1). Its use allows to 
assess whether the NNS distribution of a series of lev¬ 
els agrees with the prediction of RMT assuming that all 
levels have been observed. It is given by the expression 

/^Brody(s, /3) = , (I) 

where a(/3) = T ■ Reference [T7] reported a global 
study, involving all data points at fields B > 30 G, which 
pointed to an intermediate behavior between Poisson and 
WD. To analyze the chaoticity of the system as a func¬ 
tion of the external magnetic field (or excitation energy), 
and in particular the possibility that chaos sets in only 
starting from a certain value of the magnetic field, we 
have also performed an energy-resolved study: We have 
divided the spectra in windows of 60 levels around a par¬ 
ticular magnetic field and calculated the corresponding 
Brody parameter. We note that sixty levels is not a 
large enough number to retrieve strong statistical evi¬ 
dence from the fit, as we illustrate by means of numer¬ 
ical experiments in the Appendix!^ Nevertheless, cal¬ 
culation of this moving fit to Eq. (W provides us with a 
closer look at the behavior of the system that will permit 
a more complete understanding of the results that follow. 

The results of this analysis of short-range fluctuations 
are shown in Fig.[^ For the case of ^®®Er, we see that /? « 
0.5 throughout the data set. This seems to point that 
this isotope is in an intermediate case, not featuring full 
chaos, although it is difficult to arrive to any conclusions 
taking into account the errors of the fit. The case of 
lesEr seems to be clearer. For magnetic fields R < 20 G 
the system has not reached full chaos and presents an 
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intermediate behavior between chaos and regularity (/3 « 
0.7). Then, for i? « 30 G, we get ^ « 1 and the system 
appears to reach full chaos. However, for larger values 
of the magnetic field, the results are hindered by some 
different problem as /3 falls back to « 0.5. It is hard to 
imagine any physical mechanism inducing a decrease of 
chaoticity as the magnetic field is increased. We think 
that these results are consistent with some small loss of 
levels for high values of the magnetic field. Indeed, it is 
known that an incomplete data set from a fully chaotic 
series will show intermediate features. The analysis of 
the long-range correlations shown in the next subsection 
is consistent with this preliminary conclusion and allows 
us to be more quantitative. 

C. Long range correlations 

In order to make a complete analysis of the chaotic¬ 
ity of a many-body system, long-range correlations are 
as important as short-range correlations [7]. There are 
many different statistics to study long range correlations, 
they are usually relatively complex to calculate and dif¬ 
ficult to analyze statistically. A very useful one because 
of its simplicity is the power spectrum of the Sq statistic. 
This is defined as 

= ~ = ^ 9+1 “ ’ ( 2 ) 

i=l 

from which the power spectrum is introduced as [251126j 



where L is the total number of levels in the sequence. As 
shown in [25l|27], Ps(k) oc 1/k for chaotic systems while 
Ps{k) oc 1/fc^ for regular systems, which allows a clear 
distinction between these two effects. 

The results of our numerical studies of the long-range 
statistics on the experimental data are shown in Fig.j^as 
symbols. In Ref. [20], a very complete study of the effects 
of missing levels and mixed symmetries on the long-range 
correlations of the spectra was made and general formu¬ 
las for Ps{k) were derived. Under some simplifying as¬ 
sumptions, these formulas can be adapted to fit relatively 
small spectra (of the order of 100 levels), as is the case 
with the present erbium data. These assumptions are 
that the different independent symmetries mixed in the 
spectra have roughly the same number of levels, that the 
missing levels are distributed uniformly along the spectra 
and for the different symmetries, and that the statistical 
spectral properties of each independent symmetry with¬ 
out any levels missing coincide with the results of the 
GOE (i.e., each series is fully chaotic). Then, one has a 
two-parameter formula to fit [see Eq. (28) in |20|]. The 
two parameters are I, the number of mixed symmetries, 
and (j), the observed fraction of levels. As discussed in 


[20], as long as the original sequences follow GOE statis¬ 
tics, the values for these parameters as obtained from the 
fit are good estimations of the real parameters even if the 
other assumptions of the formula are not strictly fulfilled. 

Figure [^ shows the results of this analysis as lines. 
It follows that the data for ^®®Er at fields H > 30 G, 
and for ^®®Er at H > 25 G, are consistent with having 
only one symmetry present in the spectrum and a frac¬ 
tion of observed levels of (^ = 0.83 ± 0.08 for ^®®Er and 
(j) = 0.75 ± 0.08 for ^®®Er. The formula for the power 
spectrum as a function of (j) in the case of only one sym¬ 
metry present is written in the Appendix, Eq. ( ]AI] ). The 
presence of only one symmetry is expected, as collisions 
between erbium atoms occur at sufficiently low tempera¬ 
tures (T < 1 mK) that only the s-wave channel is open. 
At the same time, due to the substantial anisotropy of 
the interaction potentials between lanthanide species [21- 
I30j . collisions couple many possible internal states, with 
large rotational excitations, that need to be taken into 
account leading to the complex spectra observed. 

In summary, the statistical analysis of the level posi¬ 
tions point to ^®®Er featuring a chaotic spectrum at fields 
H « 30 G, while an apparent 20 — 25% fraction of missing 
levels precludes further conclusions at larger fields or for 
^®®Er. In the following section, we will perform an inde¬ 
pendent study of the level widths of both isotopes, and 
determine whether this enables to estimate the missing 
fraction. 

In principle, it is possible to estimate the fraction 
of missing levels from the short-range fluctuations as 
well m- Actually, with high-quality data, it should be 
possible to distinguish between a P{s) intermediate be¬ 
tween GOE and Poisson due to missing levels and inter¬ 
mediate statistics with a better fit to a Brody distribu¬ 
tion. Yet in practice, as we illustrate in the Appendix, 
short-range fluctuations are very ill-suited to fit a frac¬ 
tion of missing levels for small level sequences while long- 
range fluctuations, including the power spectrum, work 
very well even with sequences as small as 100 levels. 


III. STATISTICS OF RESONANCE WIDTHS 

We proceed in this section to perform a statistical anal¬ 
ysis of the observed resonance widths. It is one of the 
tenets of quantum chaos theory that “the widths of suc¬ 
cessive resonances are statistically uncorrelated” m , in 
other words, the positions and widths of levels are in¬ 
dependent quantities, containing qualitatively different 
information of the physical system [3]. For the case of 
a fully chaotic system, partial widths. A, for the decay 
through z/ channels are expected to be distributed ac¬ 
cording to a distribution with v degrees of freedom. 
This implies that level widths for a series of levels with 
the same symmetry (v = 1) follow the Porter-Thomas 
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(top) and ^®®Er (bottom). Symbols stand for data from ex¬ 
periments, while lines correspond to the Poisson prediction, a 
fit to Eq. §, and the GOE prediction (from top to bottom 
in each panel). 

(PT) distribution [5^ . 

where A is the average width. This prediction has been 
confirmed in a variety of experiments in nuclear |33j . 
atomic [31], and molecular spectra [M [551, as well as 
with microwave cavities |3il37|. 

Of particular interest for our present study is the pos¬ 
sibility to estimate the fraction of missing levels from 
a comparison of the distribution of level widths with 
Eq. ([41) [11131]. Indeed, if only widths of magnitude 
A > Amin are measured, one expects to have an observed 
fraction given by 

0PT = / PpT(A)dA = erfc ■ (5) 

The dependence of this estimate on A^iin is shown in 
the upper panel of Fig [^ Here, we have also included 
an estimate of the observed fraction, using the value of 
the smallest width reported in m to estimate Amin = 
4.3 mG, and the average of observed widths as an esti¬ 
mate of A = 30 mG [33]. From these values, it follows a 



^min/^ 



FIG. 4. (Golor online) (top) Observed fraction estimated 
from the PT distribution Q if only widths A > Amin are ob¬ 
served. The full line is the analytical estimate the symbol 
stands for the estimate of the situation corresponding to the 
experiments in Ref. |17| (see text for details), (bottom) Prob¬ 
ability density of widths, T’(A), for ^®®Er. Probabilities de¬ 
rived from experimental data are shown as shaded bins while 
lines correspond to best hts to the PT Q (dashed line, x?ed = 
5.634) and MPT <§ (solid line, Xied = 3.116) distributions. 
The hashed rectangle illustrates the range of minimum exper¬ 
imental magnetic field resolution, SBres = 5—10 mG [IT]- The 
inset shows analogous results for ^®®Er with Xred(PT) = 6.842, 
X?,d(MPT) = 3.572. 


rough estimate i^pt = 71%, which is in reasonable agree¬ 
ment with the value obtained in Sec. m based on level 
positions. 

In the bottom panel of Fig. [^ we show fits of the dis¬ 
tribution of experimental widths of both Er isotopes to 
Eq. (ji]) as dashed lines. Due to the reduced number of 
small-width data points, the fit to the PT distribution 
has been done taking into account only widths larger 
than 15 mG. Still, there is a noticeable discrepancy be¬ 
tween the best fit to the PT distribution and the experi¬ 
mental data (as goodness-of-fit measures, values for Xred, 
the chi-square statistic divided by the number of degrees 
of freedom, are reported in the caption). Among the var- 
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ious differences, we remark: 

(i) The data do not seem to drop off as quickly as 
Eq. Q predicts. As large level widths should be 
rather scarce, this may be due to statistical fluctu¬ 
ations, that happen to feature “too many” broad 
levels in the relatively small experimental data set. 
This scenario appears to be rather unlikely as the 
calculated value of the statistic depends mainly 
on the deviations for large-width levels and the 
corresponding p values (the probability of obtain¬ 
ing the observed data from a sample following the 
fitted distribution) are smaller than 0.002 for all 
the cases, including the fits to the modified Porter- 
Thomas distribution explained below. 

(ii) It is possible that some levels with small or interme¬ 
diate values of A are missing, leading to an under¬ 
estimation of T’(A) away from the maximum ob¬ 
served widths. 

(iii) Compared with analogous results in other 
quantum-chaotic systems uniisiiiMi , we notice in 
Fig.g] a striking lack of the divergence signaling 
that the limiting behavior P(A) » const/VA 
(A <C A) has been reached. To recover this 
feature, it appears that one would need to measure 
down to much smaller widths than so far observed. 

To assess which of these effects, if any, is playing the 
most relevant role, further experimental studies will be 
needed. At this point, however, we can already note 
that the average width estimated with the available data, 
A « 30 mG, is similar to the experimental magnetic field 
resolution {SB^es = 5 — 10 mG) [T7]. Experience from 
similar studies in other contexts points that one needs 
a width resolution of order « 10“^A to make a faithful 
comparison with Eq. ([4]) [iniEiiiss]. Hence, a study of 
point (iii) would require new experimental measurements 
with a magnetic field resolution of about 10“^A = 0.3 
mG. Indeed, we note that experiments with another lan¬ 
thanide species, Dy, already observe an increase in the 
number of resonant features when the magnetic-field res¬ 
olution employed is reduced from 5—10 mG to 1.3 
mG [30] ■ Our estimate based on Eq. ([^ also points in 
the same direction. 

An alternative method to estimate the values of Amin 
and A and, hence, the missing fraction, relies on the 
modified Porter-Thomas (MPT) distribution introduced 
by Flambaum et al. |38j in their study of electromagnetic 
El transitions in Ce atoms. It reads: 


r 0 


Pmpt(A) = j cPpt{A) 

[l (V) 


( 6 ) 

with C a normalization constant. Fits of the erbium data 
to Eq. (§ are shown as full lines in Fig. They appear 
to be a much closer representation of the data, especially 


for widths smaller than the average width, which is re¬ 
flected in the smaller values of than those of fits to 
the PT distribution. This puts the focus of missing res¬ 
onances on those of small width. However, the highly 
nonlinear character of the MPT distribution, and the 
limited statistics available, make it difficult to extract 
robust information from the fitted parameters. For ex¬ 
ample, the fit to the ^®®Er data yields Amin = 9.3 mG, 
Aflt = 1.3 mG. It is suggestive to note that the fitted 
value of the parameter Amin is close to the experimental 
magnetic field resolution 53-^^^ « 5 — 10 mG. An inte¬ 
gration similar to Eq. ([^ using Pmpt with the fitted pa¬ 
rameters, and setting as the lower limit of integration the 
largest value of the experimental resolution SB,-es, allows 
one to estimate the fraction of observed levels with widths 
A > max(5Hi.es): 4>mpt = PMPT(A)dA = 95.5%, 

i.e., less than 5% of levels with widths A > dB^^s would 
have not been detected. On the other hand, as Amin > 
min((5Hi.es) = 5 mG, according to the MPT model, if all 
data had been collected with the finer magnetic-field res¬ 
olution, then we would not expect to have missed any 
level with A > SBj-gg, and the whole 20 — 25% missing 
fraction estimated in Sec, jll G] would originate in narrow 
resonances, A < min((5i?res)- 

From this evidence, and taking into account our esti¬ 
mation in Sec. 001 a sizable missing fraction, we conclude 
that most missing levels must have widths A < 10 mG. 
This implies that a higher-resolution experimental study 
to address points (ii) and (iii) would be a more direct 
way to illuminate the role of chaos on erbium molecu¬ 
lar resonances, rather than increasing the sample size by 
exploring higher magnetic fields [point (i)]. 

Finally, we note the considerable difference between 
the average width estimated from the experimental data, 
A, and that coming from the MPT fit, Agt. It is worth 
mentioning that A is only a factor « 1/10 smaller than 
the mean level spacing (« 310 — 330 mG depending on 
the isotope (cf. Ref. (HI). If this were the case, one would 
expect to have overlapping resonances in the spectrum, 
making its analysis more complicated. A recent study 
of this problem has been reported in |41] which shows 
how the PT distribution needs to be modified in this 
situation. An analysis of the erbium data following this 
approach, however, points to a rather weak coupling {g « 
3.3) so that differences between the PT distribution and 
the Pm=i(A) distribution in HD are only expected for 
widths > 15 times larger than the average width. The 
erbium data only feature a handful of such points, so 
that both fits are dominated by resonances with A « A 
and yield essentially the same results. Because of this, 
one cannot draw any hard conclusions on the existence 
of such resonance overlaps with the available data. 


IV. CONCLUSIONS 

In order to study experimentally a quantum system 
from the point of view of chaos and spectral statistics, a 
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combination of tools is needed to address the questions 
of chaoticity versus regularity taking into account possi¬ 
ble experimental errors both in missing levels and mixing 
of symmetries. Only results that are consistent through¬ 
out different statistics allow to distinguish with confi¬ 
dence between experimental issues and different degrees 
of chaoticity. We believe that the analysis with combined 
tools performed here should become a standard approach 
for the statistical analysis of experimental spectra. 

We have performed such a comprehensive, multi¬ 
faceted analysis for the molecular resonances in ^®®Er and 
lesEr observed in collision experiments at ultracold tem¬ 
peratures We have analyzed both short range and 
long range fluctuations in the resonance positions, and 
the width distribution of the resonances. The conclu¬ 
sion from the spectral fluctuations is that there is a rapid 
transition to full chaos as a function of the magnetic field 
for ^®®Er, and the data are consistent with a fraction of 
missing levels of 20 — 25%. A study of the distribution of 
experimental widths in comparison with the predictions 
of RMT including a minimum observable width as a fit 
parameter, yields an independent estimate of about 5% 
missing levels above a minimum width about 10 mG; this 
highlights that the majority of missing levels must have 
widths below the reported experimental resolution. 

In summary, the present study consistently points to a 
sizable number of narrow resonances having gone unde¬ 
tected in the experiments, and provides specific recom¬ 
mendations for the magnetic field resolution needed to 
obtain higher-quality data, which are in line with exper¬ 
imental observations of molecular resonances in another 
rare-earth atomic species |40) . 
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Appendix A: Effect of small sample size on the 
analysis of spectral statistics 

As an illustration of the effect that a reduced sample 
size {L « 100 levels) can have on the estimation of the 
Brody parameter, we show in Figs, [^a-c) numerical re¬ 
sults of a fit to Eq. Q of three sets of L = 120 data points 


constructed by first creating a sample of 150 points drawn 
randomly with the WD distribution, from which then 
20% of points (uniformly distributed within the sample) 
are discarded. The data set on the first panel is best 
fitted with /Sgt = 0.90, i.e., it appears to be almost fully 
chaotic. The central panel estimates Pat = 0.81, while 
in the right panel Pat = 0.43. The figures also contain 
fits according to the procedure devised in Ref. [TS] , which 
allows to estimate the fraction of missing levels assuming 
the data correspond to a fully chaotic series. Using this 
method, we have estimated for each case observed frac¬ 
tions of pat = 95%, 78%, and 34%, respectively. Note 
how, especially for the last data series, the two fitting 
functions show greatly different behaviors—an effect here 
purely due to the reduced sample size. 

The lower two panels show the results of a similar 
study with larger samples, L = 3000 [Fig.j^d)] and 5000 
[Fig-i e)] levels, obtained by discarding 50% and 20% of 
data points from a similar initial population. These sim¬ 
ulations show that one needs to go to rather large data 
sets to be able to retrieve with confidence the missing 
fraction and/or degree of chaoticity from an analysis of 
these statistical measures. 

Finally, we performed a similar exercise for long-range 
correlations and the power spectrum. We calculated a 
number of examples with L = 135 levels coming from 
GOE matrices from which a fraction of 20% of the levels 
were taken out randomly. The power spectrum of the 
Sq statistic was calculated and results were fitted to the 
same formula as we used for fitting the power spectrum 
of the Er isotopes. 



(Al) 


where K{t) = 2t — tln(l + 2t) is the spectral form fac¬ 
tor coming from the GOE, n = L — 1 is the number 
of spacings in the sequence, and the parameter p is the 
fraction of observed levels [5D] . The results for the fitting 
parameter were always within 15% of the actual fraction 
of observed levels </> = 0.8 and for most cases the real 
fraction was within the error bars of the fit. A typical 
example is shown in Fig. 

This numerical exercise shows once more the difficulty 
to disentangle two situations (experimental missing lev¬ 
els vs. system with intermediate chaos), and highlights 
the need of high-quality data to reach reliable conclu¬ 
sions from finite data series. A more detailed quantita¬ 
tive analysis of the reliability of different statistics for 
fitting the fraction of observed levels is outside the scope 
of this work and will be published elsewhere. 
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